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1.  Introduction 


Factorial  experiments  in  which  the  additive  model  applies  have  received 
special  attention  in  the  applied  literature  because  of  the  simplicity  of 
their  analysis  and  subsequent  interpretation.  However  the  recent  emphasis 
in  the  life  sciences  on  synergism  and  antagonism  (see  e.g.  Rothman  (1974)) 
have  served  to  highlight  the  importance  of  experiments  in  which  interactions 
are  not  merely  nuisance  parameters  but  the  quantities  of  primary  interest. 

This  point  of  view  suggests  a number  of  interesting  problems.  For  example. 
Meymann  (1977)  describes  a class  of  multiple  comparison  problems  for 
determining  synergistic  effects.  This  paper  pursues  a slightly  different 
approach  and  considers  the  problem  of  selecting  the  treatment  combination 
in  a two  factor  experiment  with  both  factors  qualitative  for  which  the 
corresponding  population  interaction  is  a maximum.  The  problem  was  introduced 
by  Bechhofer  et  al.  (1977);  they  give  a detailed  analvsis  of  the  2xc  case  (two 
levels  of  the  first  factor  and  c levels  of  the  second  factor). 

The  present  work  expands  their  study  by  considering  the  design  problem  for 
arbitrary  r*c  experiments  when  the  "natural"  procedure  based  on  the  sample 
interactions  is  used. 

In  Section  2 the  model  is  introduced  and  the  design  requirement  is 
stated  for  the  problem.  In  Sections  3 and  4 the  infimum  of  the  probability 
of  correct  selection  is  studied  by  reducing  the  problem  to  a nonlinear 
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2.  Model  and  Design  Requirement 

A two- factor  experiment,  with  both  factors  qualitative,  is  to  be 

performed  with  the  first  factor  being  studied  at  r levels  and  the 

second  factor  at  c levels.  The  usual  fixed-effects  linear  model  is 

assumed  so  that  observations  Y.  (l<i<r,  l<j<c,  l<k<n) 

are  independent  and  normallv  distributed  with  means  E[Y...  ] = u..  = 

i]k  lj 

u + a . + 6 . + y . . where  T a . = F 8.  : I y..  = J y.,  = 0 are  the  usual 

i ] 11  . i H 1 .i]  - i] 

J J ! j J i J j J 

2 2 

identifiability  constraints  and  Var(Y.  ) = j"  < jo.  Here  j~  is  assumed 

i]k 

known  but  the  are  unknown.  This  paper  considers  the  problem  of 

designing  an  experiment  to  select  the  largest  (algebraic)  interaction.  Let 


7 Cl]  1 •••  1 YCrc] 


denote  the  ordered  values  of  the  (y^.K  It  is  assumed  that  the  experimenter 

has  no  prior  knowledge  of  the  pairing  of  the  Yq^ji  1 l rc  and  the 

Y . . . The  goal  is  to  select  the  treatment  combination  associated  with 
LZ 

[rc] ' 

will  be  employed  and  attention  will  be  restricted  to  procedures  which 
satisfy  the  following  probability  (design)  requirement: 


(2.1) 


P [CS]  >_  P* 


whenever  yr  , > A*  and  yr  i 
[rc]  - [rc] 


Y[rc  where  the  event  [CS] 


occurs  iff  the  treatment  combination  corresponding  to  Y^rcj  is  selected 

and  the  constants  5*,  A*  and  P*  satisfy  0 < 6*,  , ^Tyy^vr'y  K 

J (r-l)(c-l)-l 

and  — < P*  < 1 . 
rc 
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Intuitively  (2.1)  requires  that  a correct  selection  occur  with  high 

probability  only  when  both  ta>  yr  , and  ir  , , are  sufficiently 

LrcJ  [rc-ij 

far  tnart  and  (b)  yr  , is  sufficiently  positive.  When  > , , = 0 

[re]  v [rc] 

the  additive  model  holds.  The  exact  choice  of  5*,  d1'* , P*  depends  on 

economic  considerations  and  is  not  discussed  here.  The  restrictions  on 

the  and  d,Y  values  are  required  for  there  to  exist  {y^}  matrices 

which  satisfy  the  conditions  of  (.2.13  i.e.  for  the  problem  to  be  nonempty. 

The  following  "natural”  procedure,  P,  based  on  the  BLUE's,  y.^« 

of  the  will  oe  employed, 

ij 

P:  Take  n observations  on  each  treatment  combination  and 

compute  y..  - Y..  - Y.  - Y . + Y (1  < i < r,  1 < i < c) 

where  a dot  replacing  a subscript  indicates  an  average  has 

been  computed  over  the  elements  for  that  subscript.  Select 

the  treatment  combination  producing 

Yr  , i maxiY..:  l<i<r,l<j<c}  as  the  one 
LrcJ  l]  - — — — 


associated  with  y 


[rc]* 


Given  2 < r < c,  3 c and  ( 6*  ,Aa  ,Pa)  , our  problem  is  to  find 
the  smallest  value  of  n which  will  guarantee  (2.1)  when  P is  used. 

Remark  2.1.  While  procedure  P is  intuitively  appealing, its  optimality 
properties  are  unknown.  In  particular  it  is  interesting  to  note  that  a 
most  economical  property  for  P and  its  other  decision  theoretic  properties 
cannot  be  determined  from  the  works  of  Hall  1 1*3 5 9 ) , Eaton  il967)  and 
Lehmann  (1966)  since  the  joint  distribution  of  the  y.j  is  singular. 


4 


3.  Infimum  of  the  Probability  of  a Correct  Selection 

As  the  first  st-.  > in  determining  the  minimum  sample  size  to  achieve 
(2.1)  an  expression  will  be  derived  for  the  probability  of  a correct 
selection  (PCS)  when  using  P for  arbitrary  true  u. 


? [C5 1 P]  = ?[y,,  > Y . V (a,b)  * (l,i)] 

U IX  3LC i 


(3.1) 


= P 


l l Yh  • U > o,  2<_a<_r,  2 <b  < 
’ “ 2 * ;■  2 * J 


l l Y.  . ♦ l Y • >,  o,  2 < a r; 

i=2  j=2  11  j*2 

r c r 

l J *ij  + l '^ib  - 2 i b 1 c 

i=2  j*2  i=2 


since  the  sum  to  zero  over  rows  and  columns 


(3.2) 

= ?[$ 

+ y £ C]  Hf  (f),  say,  where 
^ n a. 

(3.3) 

X.  . = 
lj 

*ij  ’ Yij  (i  l1  1 i i i c)- 

(3.4) 

x’  = 

(X22 X2c,X32’,,,,Xrc)’ 

(3.5) 

Y*  = 
% 

(y«0»* • • .Y  ) and 
zz  re 

C is  a convex  cone  defined  by  (3.1)  of  the  form  {w  e R( r~i )( C~D  | Aw  oh  £ 
has  a nonsingular  (r-l)(c-l)  variate  normal  distribution  with  mean 
vector  zero  and  covariance  matrix  £ given  by 
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(3.6) 


r 1 r 2 

where  > = io..)  and  ) = (a..)  are  (c-l)*(c-l)  matrices  with 

“1  l]  L2  n 


i] 


entries 


(3.7) 


1 


(r-l)(c-l)/rcn  i = j 
-(r-l)/rcn  i i j 


and 


(3.3) 


2 

a . . 

13 


-(c-l)/rcn  i s j 
1/rcn  i i j 


and  each  row  and  column  of  \ contains  (r-1)  blocks  of 

So  the  PCS  depends  on  p only  through  ^ and  can  be  written  as 


(3.9) 


f (Y)  = PCX  £ C-y] 
no,  x x 


* / K exp{-  i £'£‘1x}dx 

C-Y 


where  K * (2ir)‘(r"i)(c'1)/2[det( £)]’1/2. 


Remark  3.1.  It  can  be  easily  checked  from  (3.1)  and  (3.2)  that  f (y) 


is  constant  under  row  and/or  column  permutations  of  the  matrix 


D 


(3.10) 


formed  from  the  vector  ^ . 

The  design  requirement  (2.1)  will  be  satisfied  if  n is  the  sma 
0 0 0 

integer  satisfying  f ( v' ) > P*  where  y = y'(n)  is  chosen  so  that 

n — % % 


(3.11) 


f (y  ) = min  f (y)  and 

n " yeF  n v 
% 


F = F(r-l,c-l,3*,A*)  is  the  convex  polytope  given  by 


*22- 


'2c 


R 


(r-i)x(c-i ) 


X ...  . .X 

r2  rc/ 


l l *„  I A*i 

i 1 


v y 


i 3 


) x . . - x , > <5 : , 2<a<r,  2<b<c; 
I ij  ab  — — 


r £ x. . + £ x . > 5*,  2 < a < r; 
i j ^ j a3  - - - 


!|vi hb a 5,>-  2 1 b s 


In  the  above  notation  row  and  column  labels  (a  and  b)  run  from  2 
and  2 to  c respectively  thus  preserving  the  analogy  with  (3.10). 
indices  i and  j range  from  2 to  r and  2 to  c respectively 


Y 


iiest 


to 

The 


A solution 


to  the  nonlinear  programming  problem  (3.11)  is  called  the 
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least  favorable  conf iguration  (LFO  since  it  represents  the  worst  case 

the  experimenter  must  design  against. 

Since  the  multivariate  normal  density  is  log  concave  and  the 

expression  (3.9)  for  the  PCS  only  involves  y as  a location  parameter 

in  the  domain  of  integration,  it  follows  directly  from  the  preservation 

theorem  of  C.  3oreIl  (1373)  (see  also  Rinott  (1976))  that  f (y)  is  also 

n -v, 

log  concave  in  y.  Hence  for  y.  ,y ^ e F and  0 < a < 1 

-v  i,-  x*.  — — 

f (&Y  + > (f  (y  )}a{f  (y„)} 

n <-*,1  — n ui  n x*. 

must  occur  at  an  extreme  Doint  of  F. 


I-a 


> min{f  (y.),f  ( y „))  and  so  the  LF 
n n -\j*. 


0 


-* . Extreme  Point  Analysis 

First  note  that  F(r-l,c-l,6*,Aft)  is  given  by 


X € R 


(r-l)«(c-l). 


T y > Aft/6*; 

ij  lj  ' 


^ y x . . / 6 * - X '6*  >1,  2 < a<  r,  2<  b<  o ; 

• ^ ab-  --  -- 


y y x.,/5ft  + y x ,/6*  > i,  2 < a < 

— — — 


1 : 


l l *H/«*  * - *ib/5*  i1-  2 1 b I c 


1 3 


and  hence  Rr-l,c-l,6ft,A*)  = 5*.F(r-l,c-l,l,Aft/6*)  where  if  a is  a 

scalar  and  S c E°  then  arS  = {ax|x  £ S).  So  it  suffices  to  determine 

the  extreme  points  of  F(r-l,c-l,l,A*)  in  order  to  solve  the  extreme 

point  problem  for  general  F since  x*  is  an  extreme  point  of 

F(r-l,c-l,5*,Aft)  x*/5*  is  an  extreme  point  of  F(r-l,c-l,l,Aft/6*). 

Next  note  that  corresponding  to  each  extreme  point  x*  of 

F(r-l,c-l,l,A*)  there  are  ( r- 1 ) ! ( c— 1 ) ! row/column  permutation  versions 

of  x*  ail  of  which  are  also  extreme  points  of  F(r-l,c-l ,l,Aft)  since 

F(r-l,c-l,l,Aft)  is  row/column  permutation  invariant.  Now  from  Remark  3.1, 

f (y)  is  constant  over  each  row/column  permutation  version  of  y and 
n \,  ‘ \, 

hence  the  LFC  y^  must  be  one  of  the  points  of  the  subset 

'U 

F'  = F' (r-l,c-l,l,A)  of  F given  by 


(4.1) 


[x  € F|x:2  > x23  > . 


- 


and 


x.„  > x,„  > 
22  — 32  — 


• -Xr2} 


r _(r-l)x(c-l) I-  . . > 

= (x  c R |Bx  > b>,  say. 
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Here  3 and  k are  of  dimensions  (rc+r+c-4)x(r-l)(c-l)  and  (rc+r+c-4)x; 
respectively. 

One  final  simplification  of  the  problem  can  be  made.  Choose  P 
orthogonal  so  that  P^P'  is  diagonal  with  entries  the  eigenvalues  of  £, 
say  A^,...,X  Then  after  the  change  of  variables  y = Px,  f ( 

becomes 


(4.2) 


f ( Y ) 
n x 


/.  . ./  K exp 
C'-Py 


, (r-l)(c-l) 

I V 


= {w  s R1"  ^]a?'w  > 0}.  Problem  (3.11)  reduces  to  determining 


Y so  that 
\ 


(4.3) 


f ( y'  ) = min  f ( y ) 
n % ycF'  n 


tor  moderate  r and  c there  are  several  efficient  algorithms 
available  for  computing  the  entire  set  of  extreme  points  of  F'  such  as 
Balinski  (1961)  and  Chernikova  (1965).  These  algorithms  systematically 
select  a subsystem  of  the  defining  inequality  system  Bx  ^_b,  solve  it 
in  equality  form,  check  the  resulting  solution  for  feasibility  in  the 
original  system  and  then  pivot  to  a new  subsystem.  There  is  an  alternative 
method  for  determining  the  extreme  points  of  F'  based  on  some  ideas 
from  graph  theory.  It  will  be  illustrated  in  the  next  section  where  the 
complete  solution  to  the  extreme  point  problem  is  determined  for  the  3x4 
case  for  arbitrary  6*  and  A*. 

While  the  enumeration  of  the  extreme  points  of  F'  is  always 
computationally  feasible,  the  same  cannot  be  said  about  evaluating  the 


I 
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objective  function  f^ty).  The  expression  (4.2)  for  f^Cy)  requires 
an  (r-l)(c-l)  dimensional  integration  over  the  convex  polytope,  C'-Py. 

In  the  3*4  example  of  Section  5,  f (y)  is  approximated  by  a simulation 

n '\j 

technique  based  on  its  probability  interpretation  (3.1).  In  summary, 

until  quadrature  or  other  techniques  are  developed  for  accurately  evaluating 

high  dimensional  integrals  over  arbitrary  convex  domains,  the  limiting 

factor  in  the  implementation  of  the  results  of  this  paper  will  be  the 

computabilitv  of  f (y). 

n -v 


w 
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5 . An  Example 

In  Bechhofer  et  al.  (1977)  a detailed  study  is  made  of  the  design  problem 
for  2*c  experiments;  closed  form  expressions  are  obtained  for  the  extreme 
points  of  F '( l,c-i ,6*, A*)  for  arbitrary  6*,  A*  and  c.  In  addition  the 
3*3  case  was  analyzed.  Tables  of  the  sample  sizes  required  to  implement 
P fo*'  .'*3,  3*4  and  3*3  experiments  are  contained  in  Bechhofer  and  Santner 
(1978). 

The  design  problem  will  now  be  studied  for  the  3*4  experiment  for 
arbitrary  6*  > 0 and  A*  > 1.26*.  First  the  extreme  points  of  F'(A*)  5 
F'(2,3,1,A*)  over  A*  * 1.2  must  be  determined.  Since  it  is  desired  to 
leave  A*  arbitrary,  the  algorithms  mentioned  in  Section  4 cannot  be  used 
to  solve  the  problem.  Instead  the  following  two-pronged  strategy  will  be 
used.  First  the  set  of  extreme  points  of  F'(A*>  above  the  face  V T x.  - A* 

i j 

will  be  determined  and  then  the  set  of  extreme  points  of  F'(A*)  on  the  face 

^ ^ x.,  = A*  will  be  determined. 

L “ lj 

The  extreme  points  of  F'ia*)  satisfying  V T x.,  > A*  can  be  found 

i j l) 

by  computing  the  vertices  of  F' (2 ,3 ,1,1. 2 ) . Now  F' ( 2 ,3 ,1,1.2 ) c F'(A*' 


since  A*  > 1.2  and  hence  everv  vertex  of 


F' (2, 3, 1,1. 2)  with  \ \ x..  * 

i 1 11 


will  be  a vertex  of  F'(A*)  above  the  face  ^ 1 1 A* . An  application  of 

i j lj 

3alins!<i's  algorithm  yields: 


(5.1) 


(5.2) 


with  T £ x..  * 1.5  and 


with  ^ ^ x 


1 ") 


Hence  for  1.2  < A*  <_  1.5  both  (5.1)  and  (5.2)  are  extreme  points  of  F ' ( A* ) 
while  for  1.5  < A*  <_  2 only  (5.2)  is  an  extreme  point  of  F'(A,l!)  and 
finally  for  2 < A*  < 00  neither  (5.1)  nor  (5.2)  is  an  extreme  point  of  F ' t A * ) . 


the  face  V £ x.. 

. = A* 

1 

• . 1 

1 

(x  c r2k3|V  y x.. 

1 *•  i-i 

, = A* 

“ “ i] 

is  more  easily  conducted  on  its  isomorphic  image,  T,  under  the  transformation 
$ given  by 


V 

A*-l 

A '--l 

A*-l  \ 

X23  X24  1 

1 -1 

A*-l 

A*-l 

A*-l 

\ 

A*-l  j 

X33  X34^ 

1 

V A!>-1 

A*-l 

A *- 1 

A*-l  / 

l V x.. 

- y x . 

- y x . * 

- y.x. 

L L ij 

L i2 

° i3 

L 14 

-Ix2j 

X22 

X23 

X24 

'^X3j 

X32 

X33 

X34 

seen 

that 

T is 

given  by 

^0 

X12 

X13 

X21 

4.  <C 

X23 

x;4  xtj>0  V (i,  j )c 

/ 

Oi 

X32 

X33 

*34  J 

l x = 4A*-3;  l x . = 4A*-4,  a = 2,3; 

j ] i 


l xu  * “*-J>  l xu>  • “‘-3-  » * 
1 1 


22  - 23  - 24’  22  - 32 
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T is  a translation  of  the  zero  row/column  sum  version  of  F"(A,l!).  Furthermore 
it  is  easy  to  check  that  x*  is  an  extreme  point  of  F"(Aft)  if  and  only  if 
$(x*)  is  an  extreme  point  of  T.  Now  T is  a (row/column  permutation  version 
of  the)  face  of  a transportation  polytope.  Transportation  polytopes  occur 
as  the  feasible  regions  of  transportation  problems  and  their  structure  is 
well  known.  In  particular  the  following  graph  theoretic  characterization  of 
the  extreme  points  of  T will  be  required  below  ( see  Klee  and  Witzgall  (1968)). 

A loop  of  x e T is  defined  as  a sequence  of  nonzero  entries  in  x 
such  that  (1)  the  row  and  column  indices  change  alternately  in  the  sequence 
and  (2)  the  first  and  last  elements  of  the  sequence  are  the  same  but  otherwise 
there  is  no  repetition. 


The  complete  set  of  vertices  of  F ' ( A * ) consists  of  (a)  those  point(s) 
(5.1)  and/or  (5.2)  which  are  in  F'(A*)  together  with  (b)  the  inverse 
images  under  $ of  those  points  (5.3)— (5.8)  satisfying  the  A*  condition. 
The  results  for  general  F(2,3,5*,A*)  can  easily  be  obtained  from  those 
for  F ' ( A ) . The  full  set  of  extreme  points  of  F'(A)  are  stated  below 
in  their  familiar  3*4  form  where  rows  and  columns  sum  to  zero: 

I.  1.25*  < A*  < 1.56*:  LFC  is  one  of 


W 


16 


Aft 

26*-2Aft 

-6* 

A ft- 5 ft  \ 

A*' -5* 

A*- 6 * 

or 

A1'*- 6* 

k 

A >'«-*$* 

26 A * 

/ 

-A*  / 

y 

A* 

26*-2A* 

A ft  - 6 ft 

-6ft  A 

y 

5A-2A'- 

A *-6* 

A *- 6 ft 

6 5* 

1 or 

1 

A*-5* 

A *-6* 

26 ft- 2A * 

0 / 
/ 

s 

4 **• 

Is 

3*- A* 

-Aft 

A*-3* 

0 

A *- 6 ft 

2 6 ft- 2A  •« 

\ 

5*-2A* 

A >'«-<$  * 

6* 

A 6 * 

J • 

After  making  the  change  of  variables  described  in  Section  4 it  is  easy 

to  check  that  f (y)  can  be  written  in  the  form  P[A(PX  t /12n/ay)  >0] 
n % ^ — 

where  X’  = < x22 ,X23 ,X24 ,X32 ,X33 ,X34 } ^ X'  ~ ( Y22  ,Y23 ,Y24  ,Y32  ,V33  ,Y34 5 

P is  an  orthogonal  matrix  satisfying 


/ 

I 

1 -2 
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A' 


\ 1 1 1 1 0 

11110  1 

1110  11 

I 

110  111 

\ 

10  1111 

y 

V\0  1 1 1 1 1 


For  anv  fixed  R = A*/5*  > 1.2  the  extreme  points  of  F ' ( 2 ,3 ,6*, A*) 
listed  in  I-V  can  each  be  written  in  the  form  6*-y(R)  where  v(R)  is  a 
6*1  vector  whose  value  is  independent  of  5*.  For  example  the  three  extreme 
points  of  I are:  <5*- (1,1, 1,1, -2,0),  <5*.  ( . 5 , . 5, . 5 , . 5 , . 5 ,-l)  and 

5*. (R-1,R-1,R-1,R-1,R-1,5-4R)  respectively.  Hence  for  fixed  R the  PCS 
evaluated  at  extreme  points  5*-^(R)  is  f(5*-^>(R) ) = P[A(PJ(  + /12n  6*/o^(R))  >_  0] 
which  is  a function  of  the  scalar  quantity  /n  6*/o.  A short  table  of  Jn  <S*/o 
values  was  constructed  for  various  P*  and  R pairs  (1.2  < R < 1.5)  by 
(a)  determining  the  zero  of  h(/n  5 */o)  = fn(6*.^(R))  - P*  for  each  candi- 
dates LFC  5*-v(R)  listed  in  I and  (b)  choosing  as  the  true  LFC  that  point 

\t 

having  the  largest  ►n’  5ft/o  value  associated  with  it.  That  largest  /n  <5*/a 
value  was  recorded  in  the  table  for  the  P*,  R coordinates.  The  values 
of  P[A(  PX  ♦ /12n  • \»(  R) ) > 0]  where  approximated  by  Monte  Carlo  Simula- 

v v • 

tion;  vectors  of  six  iid  N(0,1)  random  variables  where  generated  and  the 
set  of  inequalities  in  the  above  event  were  tested  yielding  a Bernoulli  trial 
with  the  above  "success"  probability.  The  zero  of  h(<^n  <S*/a)  was  evaluated 
by  a Robbins-Monro  type  of  stochastic  approximation  scheme. 

It  was  determined  that  for  fixed  /n  6*/o,  P[A(P£  t /12n  6rt/a*^(R))  >_  0] 
is  extremely  flat  in  R over  (1.2, 1.5);  furthermore  the  same  point. 
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5* *(1,1, 1,1, -2,0)  was  the  LFC  in  all  cases  studied.  Hence  the  following 
table  gives  values  of  /n  5*/o  which  are  accurate  for  the  entire  range 
1.2  < R = d*/5*  < 1.5  for  each  listed  P*  value. 


The  values  listed  have  been  rounded  off  from  3 place  computations  and  should 
be  accurate  to  two  places.  These  computations  were  carried  out  on 
Cornell  University's  IBM  360/168  computer. 


19 


■ 


1 


6.  Preferred  Population  Formulation 

This  section  describes  the  preferred  population  formulation  of 
Fabian  (1962)  for  determining  the  sample  size  to  be  used  with  procedure 
P.  This  formulation  is  a strengthening  of  the  indifference  zone 
requirement  (2.1)  (see  Fabian  (1962)  and  Panchapakesan  and  Santner 
(1977)). 

Fix  5*  and  A*  as  in  Section  2 and  call  the  population 

interaction  y„  (or  treatment  combination  (i,j))  preferred  (or  near 

ootimal)  iff  either  (a)  yr  , < A*  or  (b)  yr  , > A*  and  y..  > yr  ,-<5*. 

[rc]  [re]  - '[rc] 

The  goal  is  to  select  any  treatment  combination  having  a near  optimal 
y . . . Let  ft  be  the  space  of  all  r*c  matrices  with  row  and  column 
sums  equal  to  zero  and  P*  t (1/rc,  1).  The  probability  (design)  require- 
ment to  be  guaranteed  is 


(6.1) 


P [CS|P]  > P*  V Y e ft 
Y \ 


where  the  event  [CS|P]  occurs  iff  a preferred  treatment  combination  is 
selected. 

Let  ft^,  1 <_  l <_  rc  be  that  subset  of  ft  in  which  exactly  4 

Y^/s  are  preferred.  Some  ft^'s  may  be  empty,  however  it  can  easily  be 

shown  that  ftrc  = {y  e QlY[rc]  < and  ni  = W * ^Y[rc]  - A*; 

Y[rc ] “ Y[rc-1'<  - since  P^.CCS|P]  = 1 v y e ^rc  it  suffices  to 

compute  inf  P [CS|P]  for  i £ I = {l,...,rc-l}  in  order  to  determine 
°i  * 

inf  P CCS | P: . For  i e I define 

(3  l 


. 

; 


1 
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:.(y)  = ["ax{x,rc/>rrc] x 


\rc-itl  )+> [rc-itl] ^ 


> ^X(rc-i)^[rc-i]----*X(U+Y[i])] 


where  = y^j  - >[-.j  ^(i)  is  the  sample  interaction  having 

mean  y^.  ^ote  that  ? [CS|P]  = P[E;(y)]  y y € Q.  and  the  events 


I. (y ) are  nondecreasing  in  i for  fixed  y. 
1 ^ x 


Remark  6.1.  Requirement  lo.l)  imolies  that  PTE. (Y)l  > P*  V y e ; 

* 1 X — X 1 

which  is  (.2.1)  and  hence  (6.1)  is  a strengthening  of  the  indifference 
cone  approach. 


Theorem  6.1.  inf  P C -S| P]  = min  inf  P C E . (y  ) ] = inf  P[E.(y)]  and  hence 
Cl  1 i - I o Y 1 V n lx 


■"  -i  ft  £ v 

the  same  sample  sice  achieves  both  (2.1)  and  (6.1). 


?roor_.  Given  i,  2 i < rc,  it  follows  that  the  ordered  components  of 

y ; ft.  satisfy  yr  , > A*  and  yr  > ...  > yr  . , > yr  -6*  > 

x i Crc]  - [rc]  - - [rc-iti]  [rc]  - 

Yr__  n I •••  iYnv  3y  decreasing  Yr  to  Yi  . , < Yr  ,-£* 

Lrc-ij  — — lij  3 Lrc-itl]  [rc-itl]  — [rc] 

and  increasing  the  components  Yr, yr  to  preserve  row  and 

Li J Lrc-ij 

column  sums  equal  to  cero  it  can  be  seen  that  there  exists  Y'  = y'(y) 

x x v j.. 

satisfying  rc-itl  ^ l < rc  and  Y[4]lYfer  1 1 i < rc-i. 

This  implies  E.(y)  = E.(Y')  and  hence  V y e ft. 

i x i x XI 


PCE.(y)]  > PCE.(y')] 

lx  — IX 


> P[E.  . ( Y ' ) ] 

since  E.  , ( Y ' ) 

c E.(y') 

— 1-1  X 

1-1  X 

> inf  P[E . 

,(F)]  since  Y' 

~ Ted.  1_ 

x 1-1 

•lx  ^ 

i-1 

inf  P[E.(D]  > inf  P[E.  , (F)]*  2 < i < rc  and  the  result  follows  bv  induction. 
*i  1 ^ " a.  , 1-1  ** 

i-i 
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